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Abstract 

Raman intensities of Si quantum dots (QDs) with up to 11,489 atoms (about 
7.6 nm in diameter) for different scattering configurations are calculated. 
First, phonon modes in these QDs, including all vibration frequencies and 
vibration amplitudes, are calculated directly from the lattice dynamic matrix 
by using a microscopic valence force field model combined with the group 
theory. Then the Raman intensities of these quantum dots are calculated by 
using a bond-polarizability approximation. The size effects of the Raman in- 
tensity in these QDs are discussed in detail based on these calculations. The 
calculations are compared with the available experimental observation. We 
are expecting that our calculations can further stimulate more experimental 
measurements. 
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I. INTRODUCTION 



Semiconductor quantum dots (QDs) have attracted much research attention in recent 
years because of their importance in the fundamental understanding of physics and potential 
applications in electronic devices, information processing, and non-linear optics. The elec- 
tronic properties of QDs have been intensively studied in recent years, both theoretically and 
experimentally, and a clear understanding of much of the basic physics of the quantum con- 
finement effects of electrons in QDs has been achieved |jT]]. On the other hand, the vibration 
properties of QDs, i.e., the confinement of phonon modes in QDs, are less understood. 

So far, most of the theoretical understanding of phonon modes in QDs are based on the 
continuum dielectric models. The analytic expression of the eigenfunctions of LO phonons 
and surface optical phonons of small spherical p|-[lO| and cylindrical [Tl|] QDs are derived 
and the electron-phonon interactions are calculated. The extended continuum dielectric 
model [p|-[TO[] coupling the mechanical vibrational amplitudes and the electrostatic potential 
has made major improvements over classical dielectric models in the study of phonon modes 
in QDs. However, one of the basic assumptions of all dielectric models is that the material 
is homogeneous and isotropic, that is only valid in the long wavelength limit. When the size 
of QDs is small, in the range of a few nm, the continuum dielectric models are intrinsically 
limited. 

Many optical, transport, and thermal properties of quantum dots are related to phonon 
behavior in QDs. The theoretical treatment of these properties requires a reliable descrip- 
tion of phonon modes and electron-phonon interaction potential in QDs. One of the major 
difficulties of the microscopic modeling of phonon modes in QDs is its computational inten- 
sity. For example, for a GaAs QD of size of about 8.0 nm, there are 11,855 atoms in it. 
Considering the three dimensional motion of each atom, the dynamic matrix is in the order 
of 35,565. This is an intimidating task even with the most advanced computers. In recent 
years, we have developed a microscopic valence force filed model (VFFM) [p!2|-[T5[l to study 
phonon modes in QDs by employing the projection operators of the irreducible representa- 
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tions of the group theory to reduce the computational intensity. By employing the group 
theory, for example, the above matrix of size of 35,565 can be reduced to five matrices in five 
different representations of Ai, A2, E, Ti, and T2, with the sizes of 1592, 1368, 2960, 4335, 
and 4560 respectively. Therefore, the original problem is reduced to a problem that can 
be easily handled by most reasonable computers. This allows not only the investigation of 
phonon modes in QDs with a much larger size, but also the investigation of phonon modes in 
QDs with different symmetries. These investigations lead to many interesting physics that 
otherwise can not be revealed ||12|-|T5[1. With this model, we have studied the size effects 



of phonon modes in semiconductor quantum dots, including QDs of one material, such as 
GaAs or InAs, as well as QDs with a core of one material embedded in a shell of another 
material, such as GaAs cores embedded in AlAs shells. To further develop our theoretical 
model in investigations of properties of QDs, in this article, we have studied the size effects 
of Raman intensity in semiconductor QDs with the same model. 

It is well-known that the measurement of Raman spectra of a crystals is one of the most 
important methods for obtaining information about their lattice vibrations ||T^,[T^. Raman 



spectroscopy has been used to investigate the geometry, the nature, and the structures of 
QDs P,P|,p!B|-P^ . So far, most of the theoretical calculations on Raman scattering are based 



on phenomenological model Recently there are calculations based on microscopic 

models [^,^, but because of the limitation on the computational intensity, the size range 
of the QDs that can be handled is limited. 

In this article, we have calculated the Raman intensities of nanoscale Si QDs by using the 
results of the VFFM together with a bond polarizability approximation (BPA) p3|^8|. The 
calculated results are then compared with the available experimental data. We also hope 
that our calculations can further stimulate more experimental measurements on Raman 
intensities of semiconductor QDs. 

The article is organized as the following. In Sections II, we describe the theoretical models 
of VFFM and BPA; in Section III, we show our calculated results of Raman intensities and 
have some discussions; and Section IV is a summary. 



II. THEORETICAL APPROACHES 



A. VFFM for Phonon modes in QDs 

In general, the theoretical model, VFFM, that we used to investigate phonon modes in 
QDs can be used to study phonon modes in groups IV, III-V, and II- VI semiconductors. In 
this model, the change of the total energy due to the lattice vibration is considered as two 
parts, the change of the energy due to the short-range interactions and the change of the 
energy due to the long range Coulomb interaction: 

AE = AE, + AE^ (2.1) 

where the short-range interaction describes the covalent bonding, and the long range part 
approximates the Coulomb interactions of polar semiconductor compounds. For the short 



range part, we employed a VFFM with only two parameters as the following [29 



Ai?. = E^^o(^f + E^^i(A^.)' (2-2) 

where Cq and Ci are two parameters describing the energy change due to the bond length 
change and the bond angle change respectively. The summation runs over all the bond 
lengths and bond angles. Because each of these two parameters has a simple and clear 
physical meaning, this model allows us to treat the interaction between atoms near and at 
the surface appropriately. It can be further used to treat the effects of surface relaxations and 
reconstructions on the vibrations if necessary. The parameters Cq, Ci used in our calculations 
for Si are 49.1 and 1.07 eV respectively [^. Since silicon is a homopolar semiconductor 
that all atoms are neutral, the long range Coulomb interaction is not necessary. 

When considering the interaction between atoms, special attention is paid to atoms near 
the surfaces of the QDs. More specifically, for the short range interaction, when an atom 
is located near the surface, interaction from its nearest neighboring atom is considered only 
if that specific nearest atom is within the QD, and interaction from its second neighboring 
atom is considered only if that specific second neighbor atom is in the QD as well as the 



nearest neighboring atom that makes the hnk between them. The second point is important, 
because it makes sense with the physics meanings of these two parameters, but it is easy to 
be neglected. 

As we discussed in above section, we have employed the projection operators of the 
irreducible representations of the group theory to reduce the computational intensity [pi|-|5^ 



when calculating phonon modes directly from the dynamic matrices. When the results of 
phonon modes are used to calculate the Raman intensity of QDs, the advantage of applying 
the group theory to calculate phonon modes in different symmetries is even more obvious. 
Because of the symmetry dependence of Raman intensity, only phonon modes that are 
Raman active in that specific symmetry are necessary to be considered. This further reduces 
the amount of calculations required. 

B. Bond polarizability approximation 

With the vibrational normal modes obtained from the VFFM described, we use the 
BPA |23[ to calculate Raman scattering intensity of QDs. This BPA model was used in the 



prediction of Raman intensities of semiconductor superlattices [^J-|26|, fuUerences ||2^, and 
nanotubes . 



The BPA [23 1 associates an axially symmetric polarizability tensor with each bond as: 



P {Rij) = a{R,j) 7 +-f{Rij) [RijRij - ^ 7] (2.3) 

where / is a unit matrix, a{Rij) is the mean polarizability, and 'j{Rij) describes the 
anisotropy of the polarizability. Both a{Rij) and 'y{Rij) are functions of the bond length 
Rij and not direction dependent. Here Rij is the bond vector connecting atoms i and 
j, Rij is the length of Rij, and Rij = Rij /Rij. With phonon vibrations we have 
Rij = Tij + Uij = Tij + Uj — Ui, where Vij is the bond vector at equilibrium, Ui is the 
displacement of atom i, and Uij is the relative displacement of atoms j and i. For phonon 
vibrations, the condition that Uij <^ Vij always applies. Therefore, the polaribility tensor 
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P (Rij) can be expressed in power of the displacements Uij in a Talor's expansion. The 
constant term of this expansion, i.e., P (rij), can be ignored, since the total contribution of 
this term from all bonds to the Raman intensity vanishes. Keeping the first order of Uij in 
the Taylor's expansion, the polarizabihty tensor P (Rij) can be simplified as the following: 

P {Rij) w {uij ■ fij)[a'{rij) 7 +j'{rij){rijrij - ^ 7)] 

— * — * 

where a'{rij) and ^'{rij) are derivatives of a{Rij) and ^{Rij) with respect to uij evaluated 
at the equivalent bond- vector Vij respectively, and fij — Vij/rij. 

In the zincblende structure, there are four different bond orientations, so we only need 
to calculate P (Rij) for these four different types of bonds. We can choose the following 
four directions as the bond orientations: 



roi = ^(1,-1,1) 
ro2 = ^(-1,1,1) 
^03 = ^(1,1,-1) 
ro4 = -^(-1,-1,-1), 
and the four polarizations associated with the four bonds are 



'^10 0^ 



P {Roi) = —7^{Uqi,x - Uoi,y + Uoi^z)[a 



V3 



1 
1 



^ -11 ^ 



3 ' 



-1 -1 
1 -10 



(2.5) 
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3^ 



^U01,x + 2ri01,2/ - 2Moi,z -U01,x + U01,y + 2moi,z M01,i + 2^01,2/ + Uoi^z 
-Uoi,x + Uoi^y + 2^01,2 -2'Uoi,a; " 4^01,2/ " 2^01,^ 2Moi,a; + «01,2/ " ^Ol.z 



U01,x + 2lioi,2/ + '^01,2 2^01,0; + 1*01,2/ ~ '^01,2 

/ 



2lioi,a; + 2Woi.« + 4U, 



01,2/ 



*01,z 



'^10 0^ 



P {R02) 



(-U02,x + U02,y + Uo2,z)W 



1 

1 



-1 -1 
-1 1 
-11 
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7 



3v^ 



-4m, 



02,x 



2u 



02,y 



02,z 



2M02,x + 4M02,y " 2Mo2,2 2Mo2,x + ^02,5, + ^02,; 



U02,x + 2M02,y - U02,z 2Uo2,x + UQ2,y + Uq2,z '2Uo2,x " 2Mo2,j/ + 4Mo2,^ y 



(2.7) 



03 J 



1 
1 



3 ' 



1 -1 

1 -1 
-1 -1 
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3^3 



4M03,a' - 2M03,j/ + 2Mo3,^ M03,a' + U03,y + 2^03,^ -M03,x + 2Mo3,j/ + M03,, 

M03,x + %3,s/ + 2Uo3,2 -2'Uo3,x + 4Uo3,j/ + 2Uo3,^ 2Uo3,x - Uo3,y + Mo3,^ 



^ ~"'W03,x + 2Mo3,y + ^03,2 2'Uo3 .J. — ^03,^ + "^03,2 



(2. 



03,2 



04 J 



'^10 0^ 



1 
1 



3 ^ 



1 1 
1 1 



— 4Uo4,2: + 2Mo4,y + 2Mo4,2 — M04,a; " Mo4,y + 2Uo4,2 — M04,2: + 2Mo4,y — ^04,2 ^ 
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3^3 



—^04,0; — UOi,y + '^U04:,z '^UOi,x — 4'Uo4,j/ + 2^04,2 2-Uo4,x " Mo4,j/ " "^04,2 
^W04,a' + 2Mo4,j/ — Wo4,2 2Mo4,x- ~ Wo4,j/ " "^04,2 2Mo4,2,. + 2tto4,y — 4:U04,z 



(2.9) 



When we calculate the polarizability tensor for the QDs, we sum the polarizability tensor 
associated with each bond in the QDs, and get a scattering tensor 



(2.10) 



From the symmetry property, we know that the Raman tensor in cubic crystals takes 
the following forms p5|| : 



Pai-- 



^ a 0^ 



a 
0a 



(2.11) 



^ b ^ 



Pe-- 



6 
-2h 



^ ^ 



Pt2= 





d 



v 








(2.12) 



ci 
d 




(2.13) 



ci 



In our calculations we found that for semiconductor QDs with zincblende structure, after 
we applied the projection operators, P is a constant times a unit matrix for an Ai mode, 
a traceless diagonal matrix with two equal matrix elements for an E mode, and a traceless 
matrix with only non-diagonal matrix elements for a T2 mode. The Raman tensors P for 
the QDs can be expressed as the following forms: 



^ 





A; 



Pe- 



^ 2Ei ^ 
-El 
-El 



(2.14) 



(2.15) 



where 



'^0 ^ 








= + PL + pDm', 



(2.16) 



(2.17) 



Ef^ 



pLf + (PL - PLf + (PL - PL 



I \2i 



(2.18) 



and 

= [{Pu? + (PL)' + (PLn (2.19) 

where P-j is the element of the change of the polarizabihty resulting from vibration of the 
QDs in mode / as defined in (2.10). 

The Af, Ef, and are invariants of the Raman tensor and / is an index of the modes. 
This is consistent with the Raman tensor in cubic crystals of above (2.11 - 2.13). Here Af 
is related to the intensity of an Ai mode in polarized Raman geometry, and Ef and are 
related to intensities of E and T2 modes in depolarized Raman geometry respectively. For E 
or T2 modes that are two- or three-fold degenerated, the contributions of the degeneracies 
are considered in the total scattering intensity [|T^,|T^]. 

For unpolarized incident light with frequency Q scattered perpendicular to the direction 
of propagation, the intensities of Raman scattered components with frequencies Q ±uji are 



the following ^3 



(n, + l/2±l/2)^, . , 



g{uJi)[7Gf + 45Af] (2.20) 



\ /9 -\-^ /9\ 

9i^i)[QGf] (2.21) 



(nz + l/2±l/2) , ,.^^2 



2uJi 

where ni is the average occupation number of phonon mode /, and /|| and I± are intensities 
of scattered light with polarization parallel and perpendicular to the plane of scattering. In 
these equations, Af is the same as defined in (2.17), and 

Gf = 9Ef + 3T^. (2.22) 

The Raman scattering matrix for cubic crystals with arbitrary orientations and arbitrary 
incident and scattering light wave vector are studied by using Stokes vectors . The Stokes 
vectors are defined as the following: 

Let the incident elliptic polarization be described as 

E = {a + ib)EA + {c + id)EN, (2.23) 



where Ea and E^q are polarization vectors in and normal to the scattering plane. These two 
directions are labeled as A and N respectively. Its Stokes vector is defined as 



V 



(2.24) 



^ a2 + 62 + c2 + ^2 ^ 
+ 62 _ c2 - (p 
2{ac + bd) 
2{ad — be) 

For example, for a left handed circularly polarized light viewed along the direction of 
propagation of the light, its polarization vector is E = -^Ea + '^E^, and its equivalent 
Stokes vector is (1, 0, 0, -1). From this, the Stokes parameters of the scattering light are 
derived choosing four different combinations of values of a, 6, c, and d. This gives the 
complete scattering matrix aij. Assume that QDs orient randomly, all the possible incident 
directions should be averaged. Then the scattering matrix for Ai modes is 



i + cos2e -i + cos2e 
-i + cos2e i + cos2e 



V 
















2 cose 









2 COS e 



(2.25) 



and for both E and T2 modes it is 
/ 



V 



i3 + cos2e -i + cos2e 
-1 + cos2 e 1 + cos2 e 








2 cos e 

-10 cose 



(2.26) 



where e is the scattering angle. 

In our calculations, we have calculated A^, ^Ef, and 3T^ according to the right-angle 
scattering equations (2.20-2.22). For different incident light and scattering configurations, 
the Raman intensities are the linear combination of these three according to their Stokes 
vectors. 
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III. RESULTS AND DISCUSSION 



In this article we report calculated reduced Raman scattering intensity, /{rii + | ± |), 
for Si QDs with diameters from 15 to 76 A. Since a', r~-^7, and 7' are not accurately known 



a priori, these parameters used in our calculations satisfy |23 

3 3 

rih^'^ij) = gVl^j) = -^a{rij) (3.1) 

We noticed that in the early calculations of amorphous Si a' ~ was assigned 
because the observed intensity profiles I\\ and I±_ of Si have the same shape. If a' = 0, there 
would be no contribution to Raman intensity from Ai modes. However, when the size of 
QDs is small, the contribution from Ai modes might be important to consider, so we have 
assigned a' a number as above. The choice of this number is not critical, since for different 
choices of a', the shape of Raman intensity from Ai contribution is the same, and only its 
relevant strength to E and T2 modes is different. 

We have calculated Raman intensities of Ai, E, and T2 modes for Si QDs with approx- 
imate sizes of 15, 20, 25, 30, 35, 40, 50, 60, 70, and 76 A respectively, and the results are 
shown in Figs. 1-3. Here the Raman intensity / is calculated by the Lorentz broadening 

^ = ?(r)^ + (3-^z)^' ^^-^^ 

where Ii and ui are the scattering intensity and eigenfrequency of mode / respectively, and 
r is the half Lorentz width, which is taken as 2/ it = 0.64 cm~^ in our calculations. We 
have also listed the related important data in Table 1 that are numerically more clear. The 
data listed in Table 1 in order are the diameters of the QDs calculated (rf), the number of 
atoms in the QDs (A^), the intensity of the highest peak of Ai modes (/ai), the intensity 
of the highest peak of E modes {Ie), the intensity of the highest peak of T2 modes in the 
low frequency range (/r2(0) ^^nd in the high frequency range {I'r2{h)), the frequency of the 
first peak of Ai modes (cuai), the frequency of the first peak of E modes (uje), the frequency 
of the first peak of T2 modes in the low frequency range {uT2{i)), and the frequency of the 
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highest peak of T2 modes in the high frequency range {ujT2{h)). All the intensities listed 
above are Raman intensity per atom for easier comparison. We will discuss the important 
features of these results in detail next. 



A. Size effects of highest frequencies 

From Figs. 1-3 we see that in general, the major peaks in the high frequency range always 
have a T2 symmetry. The highest peaks correspond to T2 phonon modes with the highest 
frequencies. When the size of the QDs increases, this frequency approaches the frequency 
of the optical phonon frequency of bulk Si. Theoretically speaking, when the size of QDs 
approaches infinite, the total Raman spectrum of QDs approaches the Raman spectrum of 
bulk Si, and this will be the only peak left. 

In our previous calculations of phonon modes in QDs [p!^-[T5|, we have discussed the 



size dependance of phonon modes with different symmetries. We have learned that phonon 
modes with different symmetries have different size dependance, and Ai modes usually have 
the strongest size effect. For quantum dots of zincblend semiconductors, such as GaAs, when 
the size of QDs is big, the mode with the highest frequency has Ai symmetry. However, 
as the dot size decreases, there is a crossover of the symmetries, then the mode with the 
highest frequency has T2 symmetry. For Si quantum dots studied in the present paper, the 
mode with the highest frequency is always of Ai symmetry. However, the Raman intensity 
of Ai modes decreases when the dot size increases, so the strongest high frequency mode is 
always of T2 symmetry. 

When the size of QDs decreases, the frequency of this T2 peak decreases. To show this 
more clearly, we enlarged the high frequency range of Fig. 3 and plotted it as Fig. 4. From 
the data listed in Table 1, we know that when the diameters of Si QDs decreases from 75.79 
A to 14.11 A, the frequency of the highest Raman peak shifts from 518.3 cm~^ to 506.4 
cm~^ (the Raman peak of Bulk Si is at 518.9 cm~^ in our model). The systematic redshift of 
the longitudinal (LO) phonon peaks due to spacially confined phonon modes in nanocrystals 
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in the size range of a few nm have been observed |lT9|-pi|] , and recently it has been observed 
by resonant Raman scattering in three samples of Ge nanocrystals in the size range of 4-10 
nm |1^ . 

One more thing we notice from Fig. 3 and Fig. 4 for the high frequency peaks of T2 
modes is that not only the highest intensity peak red shifts as dot size decreases, but also 
weaker peaks appear at the same time. 

Experimentally it may be difficult to resolve all the weaker peaks because of broadening 
resulting from fluctuation in dot sizes. As a result one may observe an asymmetric broaden- 
ing of the Raman peak corresponding to the optical phonon as the dot size is reduced. This 
is indeed found in Raman intensities of Ge QDs [^. One may attempt to inteprete this as 
an indication that the quality of the dots may be poorer leading to larger inhomogeneous 
broadening as dot size gets smaller. However, from our calculations on Raman intensities of 
QDs, one can notice that the red shift of the strongest T2 Raman peak is smaller than the 
frequency spread of the weaker peaks which appear. In other words the broadening of the 
Raman peak is larger than the red shift as the dot size decreases. This indicates that the 
observed broadening in Raman measurements is not only due to the red shift of the peak 
alone, but there is also a contribution to this broadening from quantum size effect. 



B. Size effects of lowest frequencies 

Prom Figs. 1-3, we see that the frequency of the first peak in the low frequency range 
for all three different symmetries (Ai, E, and T2) increases as the size of the QDs decreases. 
From the data listed in Table 1, we see that when the size of Si QDs decreases from 75.79 
A to 14.11 A, the frequency of the first Raman peak of Ai modes shifts from 28.6 cm~^ to 
129.8 cm~^, the first Raman peak of E modes shifts from 8.6 cm~^ to 46.7 cm~^, and the 
first Raman peak with T2 symmetry shifts from 12.4 cm~^ to 58.8 cm~^. The size effects 
of lowest frequencies of phonon modes in QDs have been discussed in detail in our previous 
studies [|r^-[T5|, and again this is shown in the calculated Raman spectra. Furthermore, 
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we notice that even though the frequencies of the lowest frequency peak increase in all 
these three figures, they increase at a different rate. The lowest frequency of the Ai peak 
increases much faster than that of the other two. To show this more clearly, we plot the 
lowest frequency peaks versus the sizes of QDs in Fig. 5. It is obvious that the lowest 
frequency of the Ai peak increases much faster than that of the other two, because that the 
Ai modes have the strongest quantum confinement effects ||T^Jl3|l . 

Another feature of Fig. 5 is that the lowest frequencies of the Raman peaks in the 
acoustic range are roughly proportional to the inverse of the QDs diameters, which was first 
observed by Duval and his co-workers p|. This was recently observed in Si nanocrystals ||38|| , 



and it was noticed that the depolarized Raman spectra appear at much lower frequencies 
than the polarized ones. Not only our results agree with the experimental observations, we 
also learned that this is actually due to the symmetry dependence of the confinement effect 
of phonon modes, i.e., the Ai modes have the strongest confinement effects. 



C. Folding of the acoustic phonons 

In Fig. 1, we see that at the low frequency range of the Ai modes, the Raman spectra 
is dominated by a series of nearly evenly spaced peaks in the acoustical phonon range. 
As the size decreases, the spacing increases. This can be understood from the folding of 
acoustic phonons. Since the Ai modes vibrate in the radical direction, when the radius of 
the QDs increases approximately one lattice constant, there will be one more folding due to 
the confinement of the QDs along the radial direction. This should be observable by Raman 
scattering, and we are expecting such observations. 



D. Size effects on strength of Raman peaks 

Since Raman intensities plotted in Figs. 1-3 are in arbitrary unit, the size effects on 
strength of Raman peaks are not shown clearly in these figures. To show it more clearly, we 
have plotted the Raman intensity per atom versus diameters of QDs for the low frequency 
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Ai peaks, the low frequency E peaks, the low frequency T2 peaks, and the high frequency 
T2 peaks in Fig. 6. For Ai Raman spectra there are several high peaks, and we choose the 
strength of the highest peak (when the size is less than 30 A, this peak is not the peak with 
the lowest frequency). In Fig. 6 we see that the strength of low frequency peaks (Ai, E, and 
T2) decreases fast as the size of QDs increases, and the strength of high frequency peaks 
(T2) remains a constant in QDs of all sizes. This indicates that even though in bulk material 
only one major peak can be measured, when the size of QDs decreases, other peaks in the 
low frequency range will appear. Of these low frequency Raman peaks, the most noticable 
ones are probably the evenly spaced Ai peaks in the polarized spectra. Such evenly spaced 
Ai peaks should be observable. 

We want to emphasize that in Fig. 6, the Raman intensity shown is from one calculated 
highest Raman peak, either in the low frequency range or in the high frequency range. 
Experimentally the presence of multiple Raman peaks as shown in Figs. 1-4 may not be 
resolvable due to size fluctuation, and in stead a broadened peak is observed. Typically in 
such situation all the Raman intensity should be intergrated together to obtain the total 
strength of the Raman peak. To compare with the experimental results, we summed the 
calculated Raman intensities (without broadening) for all Ai, E, and T2 modes respactively 
and show them for QDs with different sizes in Fig. 7. We can see that all the Raman 
intensities will increase when the size of QDs get smaller. When the QDs get larger, the 
Raman intensity of T2 mode will approach to the Raman intensity of bulk crystal and others 
will approach to zero. This is in qualitative agreement with what was observed in Ge QDs 



E. Size effects on mode mixing 

One more thing we want to comment on from Figs. 1-3 and Fig. 6 is that for large size 
QDs, the major peak of the Raman spectra is derived from the T2 high frequency mode. 
This peak approaches the optical phonon peak in the bulk Raman spectra when the size 
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of QDs is large. When the size of QDs decreases, more and stronger peaks at the lower 
frequency range show up, which are derived from the Ai modes in the polarized Raman 
spectra and E and T2 modes from the depolarized Raman spectra. As can be seen from 
Figs. 1-3 and Fig. 6, the intensities of the Ai, E, and T2 modes are of nearly the same 
magnitude for small dots, which indicates the mode mixing due to the quantum confinement 
of phonon modes in small QDs. 

IV. SUMMARY 

In summary, we have calculated the Raman intensities of Si QDs with up to 11,489 
atoms (about 7.6 nm in diameter). The phonon modes are calculated directly from the 
lattice dynamic matrix of a microscopic VFFM by employing the projection operators of 
the irreducible representations. Based on the results of phonon modes, the Raman inten- 
sities are calculated by using a BPA. The size effects of the Raman intensity in QDs are 
discussed in detail based on these calculations. Our calculated results agree with the existing 
experimental observations, and we are expecting that our calculations will stimulate more 
experimental measurements of Raman intensities of QDs. 
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FIGURES 

FIG. 1. Reduced Raman intensities of Ai modes for Si QDs with approximate diameters in A 
indicated. 

FIG. 2. Reduced Raman intensities of E modes for Si QDs with approximate diameters in A 
indicated. 

FIG. 3. Reduced Raman intensities of T2 modes for Si QDs with approximate diameters in A 
indicated. 

FIG. 4. Reduced Raman intensities of T2 modes enlarged at the high frequency range for Si 
QDs with approximate diameters in A indicated. 

FIG. 5. Frequency of the lowest Raman peak of Ai, E, and T2 modes versus size of the dots 
for Si QDs. 

FIG. 6. Raman intensity per atom of the highest low-frequency peaks of Ai, E, and T2 modes 
and the highest high frequency peaks of T2 modes versus size of the dots for Si QDs. 

FIG. 7. Integrated Raman intensity per atom of Ai, E, and T2 modes versus size of the dots 
for Si QDs. 
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TABLES 

TABLE I. Raman intensities of Ai, E, and T2 modes for Si QDs with approximate sizes of 15, 
20, 25, 30, 35, 40, 50, 60, 70, and 76 A. The data Ustcd in order arc the diameters of the QDs (d), 
the number of atoms in the QDs (N), the intensity of the highest peak of Ai modes (Iai), the 
intensity of the highest peak of E modes (Ie), the intensity of the highest peak of T2 modes in the 
low frequency range (It2{1)), the intensity of the highest peak of T2 modes in the high frequency 
range {lT2{h)), the frequency of the first peak of Ai mode {loai), the frequency of the first peak 
of E mode (loe), the frequency of the first peak of T2 mode in the low frequency range {ujt2{1)), 
and the frequency of the highest peak of T2 mode in the high frequency range {uJT2{h)). All the 
intensities listed here are Raman intensity per atom. 



d 


N 


Iai 


Ie 


It2{1) 


lT2{h) 




LVE 


UJT2{1) 


UlT2{h) 


14.11 


87 


7.65 


9.70 


10.99 


9.88 


129.8 


46.7 


58.8 


506.4 


19.39 


191 


5.88 


5.58 


6.61 


9.60 


98.4 


35.2 


46.9 


511.2 


24.74 


417 


3.55 


3.53 


4.79 


9.47 


83.3 


27.2 


37.5 


514.3 


29.75 


705 


2.69 


2.76 


3.66 


9.22 


69.0 


23.3 


31.2 


515.6 


34.67 


1099 


2.46 


2.03 


2.64 


9.25 


61.4 


19.7 


26.6 


516.5 


39.91 


1707 


1.84 


1.51 


1.97 


9.24 


53.2 


16.9 


22.9 


517.1 


50.00 


3265 


1.27 


1.01 


1.34 


9.47 


43.2 


13.4 


18.8 


517.7 


59.99 


5707 


0.91 


0.72 


0.93 


9.67 


36.0 


10.9 


15.5 


518.1 


69.97 


9041 


0.68 


0.54 


0.69 


10.04 


30.9 


9.2 


13.4 


518.3 


75.79 


11489 


0.58 


0.46 


0.59 
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8.6 


12.4 


518.3 
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